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ABSTRACT 

Stars in the solar neighbourhood do not have a smooth distribution of velocities. 
Instead, the distribution of velocity components in the Galactic plane manifests a great 
deal of kinematic substructure. Here I present an analysis in action-angle variables of 
the Geneva-Copenhagen survey of ~ 14 000 nearby F & G dwarfs with distances and 
full space motions. I show that stars in the so-called "Hyades stream" have both angle 
and action variables characteristic of their having been scattered at an inner Lindblad 
resonance of a rotating disturbance potential. This discovery seems to favour spiral 
patterns as recurrent, short-lived instabilities. 

Key words: galaxies: evolution - galaxies: kinematics and dynamics - galaxies: spiral 



1 THE LOCAL DISTRIBUTION OF STELLAR 
VELOCITIES 



X 



The HIPPARCOS satellite (|ESA|1997| ') obtained accurate dis- 
tances and proper motions for ~ 110 000 stars, but did not 
determine radial velocities for anv. lDehnen fc BinnevI (| 19981 ) 
tried to construct an unbiased sub-sample for their analyses 
of local stellar kinem atics, and estim ated the missing radial 
velocities statistically iDehnenI l| 19981 ') . 

The G eneva-Copenhagen surv ey of nearby stars (here- 
after GCS, iNordstrom et ai]|2004 ) was designed to supply 
the missing velocity component for a homogeneous sample 
of 16 682 nearby F & G dwarf stars. This gargantuan effort 
has resulted in a large sample of stars in the solar neigh- 
bourhood with di stances and full space motions, the latest 
revision of which (|Holmberg et al\\200^ ) forms the basis of 
this paper. 

The greater part of the GCS sample (14139 stars) have 
multiple measurements of radial velocities to c heck for con- 
tamination by binaries. iHolmbe rg et all ( |200S| ) used the im- 
proved distances fro m the reanalysis of HIPPARCOS data by 
Ivan LeeuwenI (|2007l ). but substituted photometric distance 
estimates for those stars with trigonometric uncertainties 
> 13%. Proper m otions are from HIPPARCOS and Tycho-2 
llH0g et fflZ.ll200C ). I do not use the dispute d age estimates 



20071 : iHolmberg et al\\200'i\ ) in the present 



{e.g. iReid et al. 
paper. 

This monumental survey is uniquely valuable, since it 
is free from most of the selection biases that went into the 
full HIPPARCOS sample. Aside from a concentration of 112 
stars in the Hyades cluster, their distribution over the sky 
is remarkably uniform, with a slightly higher density in the 
declination range south of 5 = —26°. 

* E-mail: sellwood@physics.rutgers.edu 



The mac h ine readable table produced by 
IHolmberg et al\ (|2009f l includes the star positions and 
the U,V,&i W components of the star's motion relative to 
the Sun in Galactic coordinates^ Of the 16 682 stars in the 
table, 596 have no distance and 2536 have blank fields for 
the U , V &i W velocities. I also discard the 363 stars having 
distances > 200 pc from the Sun and the 112 Hyades cluster 
stars leaving 13 045 stars that I use in this analysis. 

The median distance of these selected stars is 74 pc and 
the sampl e within 40 pc of the S un is believed to be near 
complete. INordstrom et al\ (|2004l ) do not supply individual 
uncertainties for each velocity, but assert that they are be- 
lieved accurate to 1.5 km s~^, with the greatest contribution 
coming from distance uncertainties. I find some evidence to 
corroborate this claim in the analysis presented here. 

While the 112 Hyades cluster stars may seem too few 
in number to affect this analysis materially, I have omitted 
them because they are, in fact, resonant stars and would 
give a spurious boost to the significance of the main result 
(see 



iFamaev et al\ (120051 ) present a further sample of north- 
ern K and M giants with distances and full space motions, 
but rather few are close enough to have acceptably small 
velocity uncertainties. I have therefore not included their 
sample of giants in this analysis. 



1.1 Phase-space Distribution 

Fig. [T] confirms the salient features of Dehnen's (1998) re- 
sults, which were based on stars lacking the radial veloc- 



^ These Cartesian velocity components are oriented such that U 
is towards the Galactic centre, V is in the direction of Galactic 
rotation, and W is towards the north Galactic pole. 
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Figure 1. The lower panels show the distributions of Galacto-centric distance R — Rq, and the velocity components U (towards the 
Galactic centre) and V (in the direction of Galactic rotation) for all Geneva-Copenhagen sample stars within 200 pc of the Sun. The 
upper panels show three projections of the 3-D distribution; while there are no obvious cor relations of either velocity component with 
Galacto-centric distance, the distribution in velocity sp ace (top right) is highly non-uniform l|Dehnenlll99i : lNordstrom et al\\2004 ). The 
lines in the top right panel mark the LSR estimated bv lSchonrich et al\ | |201C| ). 



ity, and which had resi dual concerns ove r selection biases. 
The "startling" aspect (|Bovv et aZ.ll2009l ) is the amount of 
substructure in the {U, V) plane: the distribution function 
(DF) of nearby stars appears to be far from simple. There 
is little in the way of an underlying smooth component and 
the entire stellar distribution is broken into several streams 
||Bovv et aZ.ll2009l ). The principal streams that enter into the 
discussion below are labelled in Fig. (2] 

The features are too substantial to have simply arisen 
from groups of stars th at were born with similar kine- 
mat ics (e.g. lEggenl Il996), as confirmed in detailed stud- 
jFamaev et al\ l2007l : Isensbv et al] l2007l : IBovv fc Hog3 



120091 ). It is clear that the e ntire DF has been sculptured by 
dynamical processes. While iHelmi et al] (|2006l ) suggest that 
some of the structure could be the relics of satellite accre- 
tion, most work has focused on the effects of the bar and 
spiral arms. 

iKalnai j (|l99ll ) had previously suggested from entirely 
different data, that the OLR of the bar in the Milky 
Way might be close to th e solar c ircle. After the release 
of data from HIPPARCOS, iDehnenI (2000) also attributed 
the Hercules str e am (f eature C) to the OLR of the bar. 

Simone et al\ (|2004l ). on the other hand reproduce dis- 
tributions of stars with a similar degree of substructure by 
invoking a succession of short-lived spiral transients. Other 



mo dels tha t inclu d e both bars a nd sp irals are presented 
by iQuillenI (|2003l ). IChakrabart^ (|2007h and lAntoia et aj] 
(|2009t). 



Quillen fc MinchevI (|2005l ) associate the Hyades stream 



feature with 4-fold symmetric orbits in the potential of a 
long-lived 2-arm spiral. Such orbits arise near the ultra- 
harmonic (4:1) resonance of a finite amplitude perturbation, 
where stars make four radial oscillations in a single turn 
about the galaxy in a frame that rotates with the potential. 

In all these studies, the spiral waves are put in "by 
hand", as externally applied perturbations with proper- 
ties of their choosing. The supposed spiral waves should, 
of course, have arisen from the stars themselves as self- 
consistent collective disturbances, and therefore the prop- 
erties of the waves assumed by these authors may not be 
completely realistic. In particular, their assumed shapes and 
time-dependence may bias crucial resonant interactions with 
the stars, which are responsible for the substructure they 
seek to explain. 



2 DISTRIBUTION IN ACTION SPACE 

However, the distribution in velocity space, which these var- 
ious authors attempt to mimic, is not the best projection 
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Figure 2. Logarithmically spaced contours of the surface density 
of stars projected in the U — V-plane, as shown in Fig. [l] The 
letters indicate the approximate locations of the principal "star 
streams": A - Hyadcs, B - Sirius, C - Hercules, and D - Pleiades. 



of phase space to reveal the origin of the structure. A DF 
that is in equilibrium can be a function of the integrals only 
(Jeans theorem). The classical integrals for the in-plane mo- 
tion in an axi symmetric disc are E &l Lz, but actions {e.g. 
iBinnev fc Tre mainc 200j hereafter BT08) are an alterna- 
tive set of integrals that have a number of advantages. 

In a fully phase- mixed DF, stars of a given set of actions 
are uniformly distributed in the conjugate angle variables in 
a manner that could produce a complicated structure when 
observed just in velocity space. Even if the DF is not com- 
pletely phase mixed, so that Jeans theorem does not apply, 
its structure in action space is unaffected by phase mixing, 
and therefore should be easier to understand. Furthermore, 
resonant interactions with the spirals, or with the bar, are 
best viewed in integr al space, since the only lasting changes 
are to the integrals l|Lvnden- Bell fc Kalnai3ll97a ). In this 
work, therefore, I estimate action and angle variables from 
the observed coordinates and spatial motions of the GCS 
stars, and examine the stellar distribution in the space of 
these variables. 

For motion in a spherical potential, or the symmetry 
plane of an axisymmetric potential "1>(7?, 2), the radial and 
azimuthal actions are (BT08, p221) 



Jr = — (p Rdr, and = 



(1) 



Here R is the radial speed of the star, the integral is taken 
around a full radial period, and Lz is the specific angular 
momentum of the star about the Galactic centre. 

The angles, wr & lo^, conjugate to these actions both 
increase at uniform rates defined by wr = Qr = 2n /Tr and 
= = 2ti/T^, where the radial and azimuthal periods 
of the orbit are respectively defined in equations (3.17) & 
(3.19) of BT08. Note in the limit of nearly circular orbits, 
we have Slw, — >■ f2c, ^r ^ k, and Jr — > i/ta^, where Sic, 



K, and a are respectively the angular frequency of circular 
motion, the epicycle frequency, and the radius of the star's 
epicycle. 

I will here assume that the vertical motion of the disk 
stars is decoupled from their in-plane motion, or equivalently 
that the vertical action, J^, is adiabatically invariant. This is 
justified because the vertical oscillation frequency, v, of stars 
that do not climb to great heights above the disk plane is 
rapid compared with the motion of the star in the plane. 
This assumption also requires that there are no resonances 
between the vertical motion and rotating non-axisymmetric 
structures. For disk stars whose departures from circular 
motion are not too large, the variation in the vertical fre- 
quency around the star's orbit can be neglected, and the 
vertical action may be approximated as 



Ez 



where E^ ^ \{z'^ u'^ + v1) . 



(2) 



The approximate form for Ez (BT08, eq. 3.86) holds for 
stars that do not climb to great heights above the mid-plane. 
Roughly two thirds of the sample remain within 200 pc of 
the mid-plane, while just 5% can reach z-heights > 500 pc. 

Aside from a well-known tendency for older stars to 
have larger Jz (or Ez), I have not found the distribution 
of Jz values revealing, and most distributions discussed in 
this paper are integrated over all Jz. I mention the effect of 
eliminating stars with large Jz in t)4.3l 



2.1 Evaluation of Action-Angle Variables 

Evaluation of the in-plane action and angle variables for the 
GCS stars is impossible without knowledge of the mid-plane 
potential, <I>(ii, 0), which is poorly known for the Milky Way. 
Since the rotation curve for the Galaxy in the neighbourhood 
of the Sun may not be far from flat locally, I adopt the mid- 
plane potential 



<l>(i?,0) 



Vo In 



(i?o) ' 



(3) 



with the _Ro being the radius of the solar circle, and Vb be- 
ing the orbital speed of the local standard of rest (LSR). I 
correct for solar motion by addin g the peculiar velocity of 
the Sun (11.1, 12.24, 7.25) km s"^ ((schonrich et al\\2Qm ) to 
the tabulated ([/, V, W) components. I note below (§4) that 
the results in this paper are insensitive to this choice. 

In order to compute the four action-angle variables, 
{Jr, J<t,,WR,Wtf,), I integrate the orbit of every star in the 
GCS sample in this potential from its observed position. I 
normalize both actions by dividing them by the specific an- 
gular momentum of the LSR, Lzfi = RoVo; my adopted val- 
ues of these constants (_Ro = 8 kpc and Vo — 220 km s^^) 
determine only the scaling of the actions and other choices 
would not alter the conclusions. I choose wr = at the 
apocentre of a star's epicycle and = at the azimuthal 
location of the Sun. 

This procedure implies two major assumptions. First, 
the Milky Way rotation curve is locally flat; adjustment for 
a different potential would merely distort the distribution of 
stars in the space of these variables as points on a stretched 
rubber sheet. Second, I have assumed an axisymmetric po- 
tential, whereas the Milky Way certainly has at least mild 
non-axisymmetry: the large majority of stars in the GCS 
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Figure 3. The distribution of GCS stars in the space of the two 
actions Jji/Lz^o and J^/L^fi. The approximately parabolic lower 
boundary reflects the selection of stars from the solar vicinity. 
Aside from the general decrease in density towards large and 
a skew to lower caused largely by the asymmetric drift, the 
DF in this projection also shows significant substructure. 

sample have orbits that do not take them more than ~ 2 kpc 
from the orbit of the Sun, and most stay much closer. This 
range of distances is large enough for them to experience 
some non-axisymmetric changes in the gravitational poten- 
tial. The instantaneous values of these integrals at the solar 
azimuth must differ in a periodic manner from those for sim- 
ilar stars at other azimuths, but the ordering of stars in 
and Jr should not be affected, since the LSR also follows the 
same non-axisymmetric potential. Thus, neither assumption 
affects the topology of the distribution of stars in these vari- 
ables, and no features in the distribution could be erased or 
appear if they were corrected. Therefore, neither these as- 
sumptions, nor uncertainties in the Galactic constants, i?o 
& Vb, compromises the principal scientific finding of this 
paper. 

2.2 Results 

Fig. [3] shows the distribution of GCS stars in action space, 
excluding the small fraction (994 stars) with Jn/Lzfi > 0.05 
(i.e. with epicycle radii a > 2 kpc). The near-parabolic lower 
boundary of the distribution reflects the fact that all the 
stars in this sample are currently within 200 pc of the Sun, 
and therefore must have ever more eccentric orbits (larger 
Jr) as their guiding centre radii, or J^, differ increasingly 
from that of the Sun. 

The general decrease in density with increasing Jr is ex- 
pected in any reasonable stellar distribution, and the asym- 
metry between left and right results from the fact that the 
star density increases towards the Galactic centre; more 
stars visit the solar neighbourhood from the interior, which 
gives rise to the asymmetric drift. It is clear from this Fig- 
ure that the distribution in the space of these integrals also 
has substructure; the most notable feature is a pronounced 
overdensity of stars rising with slightly negative slope from 
J<t, ^ iz,0- 

Fig. m shows the distributions of angle variables for the 
same stars. We should expect fewer stars in the upper panel 
near wr = ±7r and near wr = 0, because such stars are 



0,4 



0,3 - 




Figure 4. The upper panel shows the distribution of angular 
phases, wr, while the lower panel shows the distribution of az- 
imuthal phases, w^, for the GCS stars. 

visiting the solar neighbourhood at respectively the peri- 
and apo-centres of their orbits. The fact that the peaks are 
not symmetrically placed near wr = ±7r/2 is interesting, 
as discussed below. The azimuthal phase distribution has a 
single strong peak near Wtt, — because the stars are all local 
to the solar neighbourhood. The width of this peak, which 
has a half- width at half-maximum of ~ ±11°, reflects the 
spread in Galactic azimuths of the guiding centres for stars 
currently passing by the Sun. 



3 EXCHANGES AT RESONANCES 

The upper panel of Fig. [5] shows the Lindblad diagram for 
a Mestel disk with an infinitesimal perturbation of constant 
pattern speed flp = 0.14, chosen arbitrarily. The full-drawn 
curve marks the locus of circular orbits in this {Lz,E) plane; 
no star can lie below this curve, but bound stars with E > 
iJcirc move on eccentric orbits in this potential. 

When the potential includes an m-fold rotational sym- 
metric perturbation with pattern speed flp, stars orbiting 
in the disc encounter wave crests at the Doppler-shifted fre- 
quency m\flp — Q^\. Resonances arise when 

m{np - n^) = iqr, (4) 

with / = at corotation, and Z = ±1 at the Lindblad res- 
onances where the guiding center respectively overtakes, or 
is overtaken by, the wave at the star's unforced radial fre- 
quencyQ The loci of these three principal resonances for ar- 
bitrarily eccentric orbits are marked by the broken curves in 

^ Note 1^1 = 1/2 at ultraharmonic resonances. 
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Figure 5. The upper panel shows the Lindblad diagram for a 
simple disk galaxy model. Circular orbits lie along the full-drawn 
curve and eccentric orbits fill the region above it. Angular mo- 
mentum and energy exchanges between a wave and stars move 
them along lines of slope Qp as shown. The dotted and dashed 
lines show the loci of resonances for an m = 2 perturbation of 
arbitrary pattern speed. The lower panel plots the same reso- 
nance lines and scattering vectors in action space, where = 
and circular orbits lie on the axis Jjj = 0. The dot-dashed lines 
in the lower panel show the Lindblad resonances for an m = 4 
disturbance having the same pattern speed. 



this Figure, which intersect the circular orbit curve where 
the more familiar epicyclic definitions of resonances apply. 
These lines are drawn for a bi-symmetric (m — 2) distur- 
bance; the Lindblad resonances would be closer to corotation 
for m > 2. 

Stars moving in a non-axisymmetric potential that ro- 
tates at a steady rate conserve neither their specific energy, 
E, nor their specific angular momentum, L^. But the com- 
bination 



Ij = E — Q,pLz 



(5) 



known as Jacobi's invariant (BT08, eq. 3.112), is conserved. 
Therefore, when a star changes its angular momentum by an 
amount AL^ caused by the rotating potential disturbance, 
it also changes its energy by the amount 



Thus changes to these classical integrals caused by interac- 
tions with a steadily rotating disturbance have slope Q.p in 

the Lindblad diagram. 

iLvnden- Bell fc Kalnai3 (| 19721 ) showed that a lasting 
change to the angular momentum of a star can occur only at 
resonances. As the slope of the circular orbit curve at corota- 
tion is Q.p, stars that exchange energy and angular momen- 
tum there do not move away from that curve, to first order, 
en abling exchanges to occur without heating, as described 
bv lSellwood fc BinnevI (|2002l ). The other vectors show that 
exchanges at the Lindblad resonances, on the other hand, 
cause stars to move onto more eccentric orbits (farther from 
the circular orbit curve) when angular momentum is redis- 
tributed outwards. This is the root cause of disk heating by 
spirals. 

The lower panel of Fig.[5]shows in action space the same 
resonance lines and scattering vectors as the upper panel 
shows for the space of the classical integrals. Exchanges at 
corotation cause horizontal shifts, implying no increase in 
Jh, while gains in at the outer Lindblad resonance, here- 
after OLR, and losses at the inner Lindblad resonance (ILR) 
both cause increases in Jh- The labelled resonance lines are 
drawn for m = 2 and an arbitrarily chosen value of Q.p; the 
diagram for other adopted pattern speeds in this scale-free 
potential would have a similar appearance, but with the J^- 
axis rescaled. The dot-dash lines in the lower panel mark the 
Lindblad resonances for an m — 4 wave of the same pattern 
speed, with the scattering vectors again shown. 

Notice that the direction of the scattering vector closely 
follows the resonant locus for m = 2 (dotted curve) only at 
the ILR. Thus, when stars are scattered at this resonance, 
they stay on resonance as they gain random energy, allowing 
very strong scatterings to occur. The opposite case arises 
at the OLR, where the star is moved off resonance by a 
small gain of angular momentum. The alignment is closest 
for m — 2 disturbances in potentials that yield fiat rotation 
curves, but scattering is always stronger at the ILR than at 
the OLR for m ^ 2 and in more general potentials. 

All three resonance lines have similar, though not iden- 
tical, slopes in the action plot and for all m ^ 2. Thus 
resonant stars will lie along a line in action space that has 
approximately the same negative slope for any resonance. 
However, we can identify the resonance at which they were 
scattered or trapped by examining the angles, since the res- 
onance will have reset the values of the phase angles for 
trapped stars to have the relation 



mw^ + Iwr = constant, 



(7) 



AE 



(6) 



where I and m were defined above. Both angles increase with 
time but the equation continues to hold at any fixed time 
because resonant stars have mw^ -f Iwr = mQ,p. Thus the 
value of the "constant" in eq. (0 varies as m^p{t — to), but 
from an unknown value at some earlier to that depends upon 
the phases of the stars relative to the perturbation. We will 
use eq. (Q in § g^l 



4 IDENTIFICATION OF RESONANCES 

While the principal feature in Fig. [3] is clear to the eye, I 
need a method to determine its statistical significance and 
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Figure 6. The solid line shows mean density of GCS stars in 
the (J0, Jij) plane along scattering trajectories from circular or- 
bits, as a function of the pattern speed of the perturbation. The 
shaded region shows the 99% confidence limits from similar size 
star samples with randomly selected Galacto-centric distance and 
velocity components. 

to search for possible weaker features. This requires a statis- 
tic that contrasts the data with a comparable featureless 
model. Ideally, a suitable featureless model would have a 
simple distribution function but the velocity distributions 
shown in Fig. [T] are not simple; the [/-distribution is not 
Gaussian (it has a large kurtosis), while the V-distribution 
is complex and skew. Thus no simple function of the inte- 
grals could yield a DF that would be an approximate match 
to the density of GCS stars in action space. 

4.1 Statistical Significance 

I therefore adopt a bootstrap approach in which I resample 
the actual distributions of 7?, f/ & V to generate compari- 
son distributions of pseudo stars. Selecting all three coordi- 
nates independently destroys the correlations that give rise 
to any features in integral space. While such a procedure 
could be dangerous if the velocity components have large 
and widely differing errors, the velocity and distance errors 
in the present sample are all small with respect to the ranges 
of the data. Samples of pseudo data selected in this way pop- 
ulate the (J^, J_R)-plane in a very similar manner to the real 
data, but with any small-scale structure smoothed away. 

As a quantitative estimator of structure, I compute 
the mean surface density of stars in (J^, Jh) space along 
resonance lines. I estimate the surface density of stars 
in the (J,^, Jj;) -p lane using an adaptive smoothing kernel 
ISilvermanl 1 19861 ) of elliptical shape with axis ratio 12:1 to 
compensate for the differing ranges of and Jr. To com- 
pute the mean density as a function of ilp, I integrate the 
estimated star density in this plane along a resonance line 
for that selected fip. The path of integration is from where 
the resonance line intersects the lower boundary of the stars 
to Jr = 0.05; dividing by the length of this line yields the 
mean density. 

Figure [6] compares the mean density along ILR lines 
for a range of Sip with 99% confidence intervals (shaded) 
of the same quantity selected from 2000 samples of scram- 
bled data. The Figure indicates one highly significant feature 
with a pattern speed Jlp ~ 0.296 ±0.005. No other feature is 
nearly as significant. Repeating this analysis for trapping at 
corotation and at the OLR, again revealed a single feature 
of high significance. 



0.4 
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Figure 7. The distribution of three linear combinations of phase 
angles as tests for correlations characteristic of resonances. 



In >i4.3l I report the results of tests for possible selection 
biases in this sample of stars that might affect the signifi- 
cance of the feature detected in Figure (6] 

Figure [S] assumed a bisymmetric spiral pattern, but it is 
possible the resonance arose from a spiral having higher ro- 
tational symmetry. Repeating the analysis assuming m — 3 
& 4 revealed maxima that were no less significant from that 
for m = 2 but having different frequencies. Thus not only 
do the action data not determine which resonance is respon- 
sible, but they do not determine the rotational symmetry of 
the pattern either. 

4.2 Phases 

We can tell which resonance was responsible for the density 
excess in the (J^, Jij)-plane by examining the distribution 
of phase angle variables for the stars. The distribution of 
wr (Fig. [4| has two peaks away from zero radians, while 
the distribution of values is narrow because all the stars 
in the sample are near the Sun. Because of this selection 
effect, we should expect the distribution of combinations of 
miO0 ± Wr to resemble that of ±wr alone, were the two 
angle variables uncorrelated. 

Fig. [7] shows the distributions, modulo 27r, of three lin- 
ear combinations of the two angle variables. Our expectation 
is borne out in the top panel, but not in the two lower panels. 
The change in shape of the distribution is quite dramatic, 
indicating significant correlations that are characteristic of 
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Figure 8. The upper panel shows the distribution of GCS stars in 
the space of the epicycle amplitude and spiral phase. The radial 
coordinate is a = (2Jji/ k)^^'^ , and the azimuthal coordinate is 
2w^ —wn. The lower panel is the same as shown in Fig.[3l except 
it uses the colours assigned in the upper panel. 

an ILR for a pattern. The lower two panels assume patterns 
having two- or three- fold rotational symmetry (m = 2 or 3). 
The peak near zero, which becomes more significant as m 
is increased, is a selection effect, while the feature of most 
interest is the second peak away from zero, which is still 
present for m = 4 (not shown). It is made up of stars having 
negative wr values near the peak shown in Fig. U that also 
have slightly negative values of w^. 

To investigate further, the upper panel of Fig. [8] shows 
the distribution of the GCS stars in polar coordinates of the 
epicycle radius a = {2Ja/ k)^^^ and polar angle 2w^ — wr 
as appropriate for stars trapped at an ILR of an m = 2 
pattern. The two peaks in the middle panel of Fig. [7] indicate 
concentrations of stars near 2w^ — wr ~ 0.15 and 1.17, or 9° 
and 67° . I have coloured the stars in these concentrations; all 
those having 2w^—wr = 67° ±7° and moderate radial action 
are coloured red, while those in the range 2w^ — wr = 9° ±9° 
are coloured blue. 



I have then used this colouring to show, in the lower 
panel, the locations of the same stars in the plot of ac- 
tions, reproduced from Fig. O The blue stars are concen- 
trated along the boundary of the lower plot, indicating that 
they are all close to the apocentres of their orbits and the 
concentration at this phase, discussed below, is probably of 
little dynamical significance. However, the red stars coincide 
with the resonance peak identified earlier. 

The boundaries of the triangle in the upper panel of 
Fig. [8] enclose 639 stars. These geometrically simple bound- 
aries are not intended to select all or only resonant stars, but 
it is likely that most of the selected stars are resonant. It is 
also likely that additional resonant stars lie outside these 
somewhat arbitrary boundaries. 

If we assume a m-arm spiral, we should change the po- 
lar angle to mw^ — wr. Polar plots for the cases m — 3 &i 
4 (not shown) are very similar to the upper panel of Fig. [8] 
except that the "red" stars are now concentrated at smaller 
angles, with almost equal significance. Furthermore, almost 
the same set of stars form the local peak of the distribution 
of mw^ — Wr for m = 3 & 4 as for m = 2, and they there- 
fore again lie in the density excess in the action plot (lower 
panel) . 

It is the correspondence between the stars of the se- 
lected phases that also have values of their actions to place 
them in the scattering peak that provides the most com- 
pelling evidence for an ILR. But once again the evidence 
cannot discriminate whether the stars were trapped by an 
m = 2 wave, or one having higher rotational symmetry. 

4.3 Further Tests 

I have conducted an extensive series of tests to check the re- 
sult in Fig. |5] The main concern is that the apparently reso- 
nant stars stand out because they share some other property 
unrelated to a resonance. 

Some 37% of the GCS stars are known or suspected 
binaries. The significance of the resonance peak is unaffected 
by including or eliminating the stars fiagged as possible or 
known binaries. 

Increasing the kernel width used to estimate the density 
in the ( J0, Ji})-plane (default = 0.02) smoothed the curve 
in Fig. |S] and lowered the maximum as expected. But it also 
reduced the ranges of the confidence intervals, so the feature 
remained highly significant. 

The principal source of uncertainty in the input data 
arises from the uncertainty in the distance to each star, 
which generally rises with distance to the star. I therefore 
repeated the analysis using stars lying within various dis- 
tances from the Sun (d < 50 pc at the most extreme) , which 
had little effect on either the significance or the position of 
the principal peak in the upper panel of Fig. [S] 

Performing the same analysis on stars separated into 
opposite hemispheres on the sky, either in equatorial or 
Galactic coordinates, did not affect the significance of the 
peak. Since the spectroscopic radial velocities could be more 
precise than those derived from the proper motion, I tried 
separating the stars into samples within 60° of the Galactic 
poles from those lying closer to the plane. The radial veloc- 
ity contributes little to the U Sz V velocity components for 
stars near the caps, but the scattering peak in this case is 
just as significant as for the remaining stars. I also tried sep- 
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Table 1. Estimated pattern speeds and corotation radii for the 
fitted disturbance. The possible systematic errors in these quan- 
tities are much greater than the small statistical uncertainties. 



Angular periodicity 






RCR 






km s"-"^ kpc"-*^ 


kpc 


m = 2 


0.296 


8.1 


27 


m = 3 


0.534 


14.7 


15 


m = 4 


0.655 


18.0 


12 



arate analyses for stars in quadrants centred on the Galactic 
centre and anti-centre, from those in quadrants centred on 
the apex and ant-apex of Galactic rotation, again finding 
little change. 

Finally, I investigated the effects of dividing the sam- 
ple using by the energy of vertical motion, Ez (eq. [5)). The 
resonance peak is most significant among the stars in the 
lowest ~ 66% of the vertical energy range, and much less 
significant in the other 33%. Such a dependence seems en- 
tirely reasonable, since stars with larger vertical oscillations 
will be much less affected by in-plane potential variations. 

These various tests have proved reassuring. The effects 
of selecting subsamples of the stars are all entirely consistent 
with physical expectations, and the likely influence of the 
small errors in the quantities tabulated in the catalogue. 



4.4 Adjustments to the LSR 

The local standard of rest (LSR), which is the speed of a par- 
ticle on a notional "circular orbit" passing through the Sun's 
position, may not even be well-defined because the Galaxy 
is at least mildly non-axisymmetric. Furthermore, the Sun's 
speed relative to the LSR is hard to determi ne, and esti- 
mate s by Binney and his col l aborators (Dohnen fc BinnevI 
1 19981 : lAumer fc BinnevI l2009l : ISchonrich et al\ |2010| ) from 



the HIPPARCOS data have changed by several km s~ at 
each revision. 

As the values of the action-angle variables depend on 
the choice of the LSR, I have experimente d with adopting 
values that differ from the latest estimate jSchonrich et all 
|2010D by ±10 and ±20 km s"^ in both U and V. Shifting 
the V-component essentially shifts the entire distribution 
in Fig. |3] to the left or right while preserving most of the 
structure. For larger revisions to the LSR, the upward ris- 
ing tongue becomes rather broader, however. Furthermore, 
the distribution of the angle variables is not changed qual- 
itatively, and 1 always obtain peaks in the distribution of 
mWif, — Wr for m = 2, 3 & 4 that are located away from the 
direction towards the Galactic centre. Thus the main find- 
ing of this paper is little affected by the precise choice of the 
LSR. 



5 DISCUSSION 

Fig. |S] provides strong evidence that a small fraction (~ 5%) 
of the GCS stars have been trapped at an ILR of a rotating 
disturbance. 



5.1 Angular Periodicity 

It is indeed surprising to find evidence for an ILR near the 
Sun, as it seems to imply a spiral pattern extending outwards 
from the Galactic radius of the Sun's orbit. However, none 
of the foregoing analysis could discriminate whether the ILR 
was that of a bisymmetric pattern or one of higher rotational 
symmetry. 

Assuming Ro — 8 kpc and Vo — 220 Im s~^, and us- 
ing the evidence that the scattering peak was caused by 
an ILR, we may determine the pattern speed of the dis- 
turbance from the frequency at which the mean density in 
Fig. [6]peaks. Table [1] gives, for each angular periodicity, the 
measured frequency of the peak, and the implied pattern 
speed and corotation radius, also assuming an exactly fiat 
rotation curve. (The statistical uncertainty in the measured 
Qp is < 2%, but systematic errors arising from the various 
assumptions are likely to be much larger.) While the Milky 
Way rotation curve may not be exactly flat, it is clear that 
corotation lies in the far outer disk for all of these patterns. 

If the disturbance responsible for the ILR is a spiral pat- 
tern (it could be some other slowly rotating disturbance), 
one would expect corotation to lie within the disk of the 
Milky Way, which probably does no t extend much be yond 
a Galacto-centric radius of 14 kpc (|Robin et aZ.lll993 . but 
see also Carraro et al. 2010). This argument would strongly 
disfavour an m = 2 spiral, but the Milky Way disk may 
be extensive enough to accommodate spirals of higher rota- 
tion al symmetry. 

iQuiUen fc MinchevI (|2005l ) calculated the appearance of 
the (U, V^)-plane in a simple galaxy model that was per- 
turbed by a large-amplitude, long-lived, bi-symmetric spi- 
ral that rotated with various pattern speeds. For a pattern 
speed low enough to place the ILR near the Sun, the phase- 
space structure they calculated did not correspond well with 
that observed, but higher pattern speeds that placed the 
Sun near the ultraharmonic resonance of their assumed spi- 
ral created a feature resembling the Hyades stream. They 
showed that the orbits responsible were 4:1 resonant peri- 
odic orbits in their adopted pattern. Unlike in their work, 
my analysis makes no assumption about the spatial form 
or lifetime of the potential perturbation, but the fact that 
they favoured 4:1 periodic orbits as the cause of the Hyades 
feature may be hinting that an m = 4 spiral perturbation is 
favoured. 



5.2 Properties of the Resonant Stars 

The upper panel of Fig. |9] shows that the 639 stars that I 
coloured red in Fig. [8] lie in the "Hyades stream" in {U,V)- 
space (feature A of Fig. [2| . Note that the highlighted stars 
do not include all the resonant stars; thus the fact that not 
all the Hyades stream stars are coloured, and that a few 
coloured stars are outside the stream may not be significant 
(see ^5.4p . It therefore seems clear that the larger part of the 
Hyades stream is caused by the ILR identified here. 

Aside from their kinematics, I have been unable to find 
any property of these 639 stars that clearly differs from the 
parent distribution of the GCS sample. They have a similar 
distribution of di stances, and in cidence of binarity, for ex- 
ample. Note that iFamaev et al\ ((2007) do find a mild, but 
significant, metallicity excess among the Hyades stars com- 



© 2010 RAS, MNRAS OOP. [111121 



Detection of a Lindblad Resonance 9 



40 



tr 

00 

> 



-40 



-80 




-80 -40 40 

Ulsr (km s"^) 




Figure 9. The upper panel reproduces the distribution of GCS 
stars in the top-right panel of Fig. [T] but with those shown in 
red being the 639 stars picked out in Fig. \E\ The lower panel 
shows the positions, in right ascension and declination, of the 
highlighted stars. 



5.3 The Age of the Resonance 

Were reliable ages for the stars available, we could use them 
to date the origin of the scattering peak, since it can contain 
stars of all ages except those that are younger than the wave 
that caused it. This is because the gas from which stars form 
is not subject to the laws of collisionless dynamics; a gas 
cloud in a resonance will experience the same gravitational 
accelerations as do the stars, but dissipative collisions will 
quickly destroy the resonant signature in phase space. 

While age estimates fo r individual stars appe ar to be 
unreliable l|Reid et aZ.ll2007l : iHolmberg et aZ.ll2007l ), we can 
at least use the Hyades cluster, which does appear to be 
resonant f i)5.2p . The cluster must have been formed before 
the perturbation that scattered it, which places an upper 
bound on the time since the resonance of 648 ±45 Myr {e.g. 
toe Gennaro et al\\200^ ). 

We can use also use the spread in the spiral phase, or 
the width of the peak in Fig. [T] to argue that the pertur- 
bation must have been recent. Assuming spirals to be tran- 
sient, the resonance will have had a width, AQp = O.lfip, 
say. Thus the width of the spiral phase will increase with 
time as A{m'w,f, + Iwr) = Q.lmQ.p{t — to). A generous es- 
timate of the width of the peak is 0.3 radians (or ~ 19°, 
over twice the width of the red triangle in Fig. [S]) and there- 
fore the time since the resonance t — to < 0.3/(0.1mf2p), 
where the inequality arises because the width can only be 
increased by measurement error. The value of fip is given by 
the position of the peak in Fig. [6l or ilp = 0.3Vo/Ro- Using 
Ro = 8 kpc and Vo — 220 km s~^ and adopting m — 2, 
I find t ~ to < 180 Myr, or the resonance occurred within 
the last Galactic rotation. This rough calculation made a 
number of assumptions; assuming a sharper resonance or 
a broader spread of spiral phases would increase the time 
scale, but the values adopted seem reasonable. It is reassur- 
ing that this estimate is substantially less than the age of 
the Hyades. I therefore conclude that the resonance is indeed 
very recent. 

Note that not only will phase mixing blur resonances 
in angle space, but we must expect resonant features to 
become less distinct in action space over time because 
stars are also scatter ed by giant molecular clouds {e.g. 
iBinnev fc Lacevlll988l '). The rate is hard to quantify , how- 
ever, since the rate of cloud s cattering is uncertain (jLacevI 
Il99ll : iHanninen fc Flvnnll200^ ) . 



pared with the mean for the GCS, although IBovv fc Hoed 
l|2009l ) find they are "barely distinguishable" from the back- 
ground population. 

The lower panel of Fig. [9] shows that the 639 red stars 
are almost uniformly distributed over the sky. Note that I 
have excluded stars in the Hyades cluster itself from all the 
analysis because it would be wrong to treat all the stars 
within a gravitational bound cluster as if they were dynami- 
cally independent. Had I included the cluster, an additional 
84 of the 112 Hyades stars would have been coloured red 
in Fig. [8l the resonant feature would have appeared yet 
more significant, and the cluster would have stood out as 
a strong concentration near {a,S) = (4'" 27™, 15°52'). Thus 
the Hyades cluster itself is resonant. 



5.4 Other "Streams" 

Since the red stars in Fig. [8] turned out to be the Hyades 
stream, I have attempted to identify in the same Figure the 
other principal features named in Fig. [S] I have coloured all 
the GCS stars lying in selected elliptical areas in velocity 
space in the top panel of Fig. 1101 and the same colours are 
used for each star in the lower two panels. Again I use red 
for the Hyades region, but here the selection deliberately 
errs on generous side to include most of the stream stars 
and more besides. 

Tracing these stars in other panels reveals that most of 
the Hyades stream stars lie in the resonance ridge in Fig. [3] 
and have the range of phases that yielded the peak in Fig. [T] 
The Sirius stream (green) appears not to be resonant, as its 
stars extend over a significant range of but not in Jr, and 
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Figure 10. The top panel assigns colours to some of the stars 
near the centres of the different streams: Hyades — red, Sirius 
— green, Hercules - dark blue, and Pleiades - cyan. The same 
colour codes are used to show the locations of these streams in 
the distribution of action and angle variables in the lower two 
panels. 



have no special concentration in phases. Both the Pleiades 
(cyan) and Hercules (blue) stars are near the low- bound- 
ary of the sample and may be parts of larger structures that 
extend outside the sample, if so their range of phases may 
be limited. This selection effect makes it impossible for this 
diagnostic to determine whether they are associated with a 
resonance. 



5.5 Implications for Spiral Structure Theory 

The identification of an ILR in the solar neighbourhood is a 
powerful discriminant between the different theories of spi- 
ral pattern generation. It has lo ng been known that spiral 
waves a re damped at an ILR (e.g. iBinney fc TremainabOOSl . 
p. 502). lBertin fc LinI (|l996l ). who argue that spiral patterns 
are long-lived waves, therefore require the ILR to be shielded 
in order to avoid fierce damping. Thus evidence for an ILR 
excludes their quasi-steady spiral hypothesis as the originat- 
ing mechanism for the disturbance that created it. While 
this evidence does not exclude the possibility of steady spi- 
rals elsewhere, it seems unattractive to invoke a different 
mechanism to account for spirals in places other than where 
we can mount a test. 

Strong ILR damping of individual patterns there- 
fore requires spirals to be continuously regenerated in 
order that they be com mon features in galaxies. While 
iToomre fc KalnaisI (|l99l') envisage sp iral patterns as recur- 
ring transient r esponses, only my idea (jSellwood fc Linll989l : 
ISellwoodllioOol ) invokes resonant scattering as the cause of 
a recurren t cycle of true instabilities. Indeed, I predicted 
(ISellwoodl '19941 that such a feature might be visible in the 
HIPPARCOS data. 



6 CONCLUSIONS 

The foregoing analysis and tests provide compelling evidence 
that the "Hyades stream" in the solar neighbourhood has 
been created through scattering at a recent inner Lindblad 
resonance. The strongest evidence is that the stars clustered 
over a narrow range of the combination of angles expected 
for an ILR are also those that lie in the highly significant 
overdensity in action space, as indicated by those coloured 
red in Fig.|8l These same stars, which are distributed all over 
the sky, make up the Hy ades stream as shown in Fig. |9l 

IBovv fc Hoel Hooi) construct a plot similar to Fig. |3l 
also using the GCS sample but restricted to the ~ 9 500 
stars within 100 pc of the Sun. Instead of using radial action 
as the ordinate, their Figure 17 uses the almost equivalent 
energy of non-circular motion, E — Ec, and they scale their 
plot to include a larger range of epicycle sizes. Because of 
the difference in scale, the feature I report here is rather 
inconspicuous in their Figure, but th eir aim was to te st for 
the feature I had claimed a hint of in ISellwoodI HooO). 

Thus only the Hyades stream (feature A of Fig. [2]) can 
be attributed to this ILR and the remaining structures in 
the local phase space distribution must be accounted for in 
other ways. It is likely that the OLR of the bar in t he Milky 
Way accounts for the Hercules stream (feature C, iDehnenI 
|2000D . but we still lack convincing models to account for the 
Pleiades and Sirius streams. Since neither feature stands out 
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in Fig. [To] (see also Fig. 17 of lBovv fc H"og3l2009l ). these kine- 
matic groupings of stars have probably not been caused by 
resonances with spirals, although they could perhaps still re- 
flect the continued presence of spirals in th is part of the disc 
with the resonances lying elsewhere {e.g. iDe Simone et al\ 
|2004| ). 

Most theoretical studies of spiral structure start with 
the assumption that the DF of disk stars is smooth and fea- 
tureless. The immediate implication of a prominent feature 
in Fig. [3] is that this assumption is incorrect. The fact that 
it was created by an ILR is strong evidence against the dis- 
turbance being a qu asi-steady spiral of the kind proposed by 
iBertin fc LinI (|l996l l. Rather, it supports the general picture 
of s piral structure generation that I have been advocating 
{e.g. ISellwood fc Linll 19891 : ISellwoodllioOol ). 

Since a spiral pattern should extend outwards from an 
ILR, it is somewhat surprising to flnd evidence for one this 
far out in the Galaxy. However, the transient spiral waves 
manifested by simulati ons of disc galax ies have a wide spread 
of pattern speeds (e.Q. ISellwoodl[l98^ 1. and do generally in- 
clude some spirals in the far outer parts of the disc. The ac- 
tual pattern speed implied by this resonance and the radii of 
corotation assuming it is a 2-, 3- or 4-armed spiral, are given 
in Table [T] the value for m = 2 is clearly lo wer than other 
estimates for spirals in the Milky Way (see iGerhardI [201ol . 
for an up-to-date summary). While the many assumptions, 
such as the slope of the local rotation curve, the value of the 
LSR, the neglect of non-axial symmetry, etc., are potential 
sources of large systematic errors in the values given in Ta- 
ble [1] I have argued that none compromise the existence of 
the resonance. If the disturbance responsible for the ILR is 
a spiral, then a multi-arm pattern seems more likely than a 
bi-symmetric one, but it may also have resulted from some 
other type of slowly-rotating disturbance in the outer galaxy, 
such as a rotating bi-symmetric distortion of the halo. 

Once the stil l more precise and m uch more extensive 
data from Gaia (|Perrvman et al\ 1200 ll ) become available, 
it will be possible to repeat this analysis with far larger 
samples of stars, and even at locations other than the solar 
neighbourhood. In particular, it should be possible to follow 
the azimuthal variation of stars in the resonance discovered 
here to determine its angular periodicity. It is to be hoped 
that other resonant scattering features will be found, and it 
may be possible to piece together the history of spiral waves 
in the Milky Way. 
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